Probing dynamics of an electron-spin ensemble via a superconducting resonator 
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We study spin relaxation and diffusion in an electron-spin ensemble of nitrogen impurities in 
diamond at low temperature (0.25-1.2 K) and polarizing magnetic field (80-300 mT). Measure- 
ments exploit mode- and temperature- dependent coupling of hyperfine-split sub-ensembles to the 
resonator. Temperature-independent spin linewidth and relaxation time suggest that spin diffusion 
limits spin relaxation. Depolarization of one sub-ensemble by resonant pumping of another indi- 
cates fast cross-relaxation compared to spin diffusion, with implications on use of sub-ensembles as 
independent quantum memories. 
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The study of spin ensembles coupled to superconduct- 
ing integrated circuits is of both technological and fun- 
damental interest. An eventual quantum computer may 
involve a hybrid architecture [1-4] combining supercon- 
ducting qubits for processing of information, solid-state 
spins for storage, and superconducting resonators for in- 
terconversion. Additionally, superconducting resonators 
allow the study of spin ensembles at low temperatures 
with ultra-low excitation powers and high spectral resolu- 
tion [5, 6]. While one spin couples to one microwave pho- 
ton with strength g/2n ~ 10 Hz, an ensemble of N spins 
collectively couples with g ens = g^fN [7, 8] , reaching the 
strong-coupling regime g ens > k, 7 at JV > 10 12 [8-10], 
where k and 7 are the circuit damping and spin dephas- 
ing rates, respectively. 

Among the solid-state spin ensembles under consid- 
eration, nitrogen defects in diamond (PI centers) [11] 
are excellent candidates for quantum information pro- 
cessing. Diamond samples can be synthesized with PI 
centers as only paramagnetic impurities. Additionally, 
samples with spin densities ranging from highly dense 
(> 200 ppm) to very dilute (< 5 ppb) are commer- 
cially available, allowing the tailoring of spin linewidth 
(7 oc A [12]) and collective strength (g cns oc VN). In 
contrast to nitrogen-vacancy centers in diamond [13] and 
rare-earth ions in Y2S1O5 [14, 15], PI centers are op- 
tically inactive, making a coupled microwave resonator 
an ideal probe for their study. However, the magnetic 
fields > 100 mT needed to polarize the ensemble at the 
few-GHz transition frequencies of circuits [16] must not 
compromise superconductivity. The freezing of all spin 
dynamics in a high-purity PI ensemble by the field would 
allow quenching spin decoherence [17] through dynamical 
decoupling [18], realizing a useful quantum memory. 
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Here, we investigate the dynamics of a PI electron-spin 
ensemble probed by controlled coupling to two modes of 
a coplanar waveguide (CPW) resonator. The resonator 
is patterned on a NbTiN film [19] withstanding applied 
magnetic fields beyond 300 mT. Three hyperfine-split 
spin sub-ensembles are clearly resolved over the temper- 
ature range 0.25-1.2 K. The collective coupling of each 




FIG. 1: (color online), (a) Schematic of the hybrid resonator- 
spin system. A single-crystal diamond piece (1.7 mm x 
1.7 mm x 1.1 mm, type-lb Sumicrystal, ~ 200 ppm N con- 
tent) is placed on top of one of four CPW resonators capac- 
itively coupled to a common feedline. The resonators are 
patterned on a NbTiN film (70 nm thick, critical temper- 
ature T c = 12.5 K) on sapphire (C-plane, 430 fim thick). 
An external magnetic field Bu is applied parallel to the film, 
along the diamond [100] direction, (b) Hyperfine interaction 
A « 94 MHz with the N host nuclear spin splits each electron- 
spin level into a triplet. Only electron-spin transitions that 
preserve nuclear spin (solid arrows) are allowed, (c) Bu tunes 
the electron-spin energy levels through resonance with the A/4 
or 3 A/4 modes of the resonator. The dashed line represents 
the thermal energy kuT/h « 5 GHz at T = 0.25 K. 
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sub-ensemble to the resonator reaches g cns — 17 MHz 
at 250 mK, in accordance with field- and temperature- 
controlled spin polarization, and extrapolates to g cns — 
23 MHz at full polarization. Using a pump-probe tech- 
nique in the dispersive regime, we measure spin linewidth 
and relaxation time. The observed temperature indepen- 
dence indicates that internal spin equilibration is dom- 
inated by spin diffusion across the mode volume [20] 
rather than spin-lattice relaxation [21]. Finally, as an 
initial test of the possible use of sub-ensembles as inde- 
pendent memories, we measure the steady-state depolar- 
ization of a sub-ensemble by resonant pumping of an- 
other. The pump-power dependence observed indicates 
fast cross-relaxation compared to spin relaxation in the 
mode volume, calling for follow-up experiments probing 
the millisecond scale. 

Our hybrid system, shown schematically in Fig. 1, con- 
sists of four resonators capacitively coupled to a com- 
mon feedline and a type-lb diamond sample placed above 
one of them. The electron-spin ensemble consists of un- 
paired electrons (spin-1/2) at substitutional nitrogen im- 
purities [Fig. 1(a)]. Each electron spin exhibits strong 
anisotropic hyperfine interaction with the host nucleus 
(spin-1). The Hamiltonian for one defect is given by 
i?N = —m a B ■ S + hS ■ A ■ I, with S and I the spin 
operators for the electron and nitrogen nucleus, respec- 
tively, mo/h — 28.0 MHz/mT, h Planck's constant, and 
A = diag(81.33, 81.33, 114.03) MHz the hyperfine inter- 
action tensor [22] [third (first, second) index parallel (nor- 
mal) to the Jahn- Teller axis] . Low-energy terms only in- 
volving I have been left out. We tune the electron-spin 
transitions with a magnetic field (-By) applied along the 
diamond [100] direction. Because all N-C bonds have 
(111) orientation and make the same angle with B«, the 
hyperfine interaction is the same for all impurities, cre- 
ating three hyperfine-split electron-spin transitions [11]. 

Measurements [23] of the feedline transmission 
|S2i|(/, near the fundamental (A/4) and the second- 
harmonic (3A/4) modes at T — 0.25 K and 1.2 K clearly 
show three avoided crossings [10], as expected for co- 
herent coupling [24] (Fig. 2). The coupling strength 
of each hyperfine transition to the 3A/4 mode is evi- 
dently stronger than to the A/4 mode, and decreases 
for both modes with increasing temperature. The hy- 
bridized dips observed when spin sub-ensembles are res- 
onant with the 3A/4 mode [Fig. 3(a)] support strong cou- 
pling (2g cns > 7, k). The absence of double dips on reso- 
nance with the A/4 mode indicate 2g ens < 7. 

We extract g ens using the model presented in Ref. 10, 
treating the spin sub-ensembles as separate harmonic os- 
cillators coupled to the resonator, but not to each other: 
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Here, A c = u — lu c is the frequency detuning between 
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FIG. 2: (color online). Transmission spectroscopy. Image 
plots of normalized feedline transmission as a function of B\\ 
and frequency near the A/4 and 3A/4 mode resonances at T = 
0.25 K (a,b) and 1.2 K (c,d). Each plot reveals three avoided 
crossings, corresponding to allowed hyperfine-split electron- 
spin transitions. Note that the frequency span in (b,d) is 10 
times larger than in (a,c). The arrow in (a) points to a flux 
jump shifting the resonator frequency. All other image plots 
shown are corrected for these rare events. 



the probe and bare resonator mode, K\ and k c are the 
resonator intrinsic and extrinsic dissipation rates, A„ = 
lj — tj mi=n is the probe detuning from the mj — n hy- 
perfine transition and 7 is the transition linewidth (as- 
sumed independent of mi). As shown in Figs. 3(a) and 
3(b), fitting the double-dip spectrum for the 3A/4 mode 
and the quality factors (Q) for A/4 mode at 0.25 K us- 
ing Eq. (1) yields collective coupling strengths g CDS /2n — 
17.0 ± 0.4 MHz and 3.9 ± 0.2 MHz, respectively [23]. 

To investigate the temperature dependence of the col- 
lective coupling strength to each mode, we measure trans- 
mission spectra at several temperatures in the range 0.25- 
1.2 K and perform the same analysis as above [14, 15, 25]. 
The results are shown in Fig. 3(c) together with the best 
fits to 



9 m s{T)=g CBa (0)JP(B h T) 



(2) 



where g C ns(0) is the zero-temperature coupling strength, 
P(B\\,T) — tanh(mo-B|| /2/cbT) the spin polarization in 
thermal equilibrium, and fee the Boltzmann constant. 
Two factors combine to make g ens (T) higher for the 3 A/4 
mode. First, P increases monotonically with the Zeeman 
energy itiqBu . Second, the bare spin-coupling strength g 
increases as ^fuJ^ owing to a larger vacuum magnetic field 
strength. The ratio 2.7 between the best-fit <7 e ns(0)/27r 
values for the 3A/4 and A/4 modes (22.7 ± 0.6 and 
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FIG. 3: (color online). Determination of the collective cou- 
pling strength g cns - (a) A vertical cut of Fig. 2(b) at By = 
272.8 mT (dashed arrow) shows Rabi-split transmission dips. 
The best fit to Eq. (1) gives g cnB = 17.0 MHz. (b) Measured 
loaded quality factor of the A/4 mode as a function of -By at 
T = 0.25 K. The best fit of Eq. (1) away from the avoided 
crossings gives g cns — 3.9 MHz. Arrows point to satellites re- 
sulting from the hyperfine coupling of the electron spin to the 
nuclear spin of 13 C atoms adjacent to some PI centers (see 
[23] for further discussion), (c) Best-fit <; cns to the A/4 (cir- 
cles) and 3A/4 (squares) modes as a function of temperature. 
Solid curves are the best fits of Eq. (2). Error bars are smaller 
than the symbol size. 

8.3±0.2 MHz, respectively) differs from the expected 
This discrepancy may be due to inhomogeneous distribu- 
tion of PI centers in the mode volume [26] (see further 
below) . 

Having characterized coherent coupling in the hybrid 
system, we now turn to using the resonator as a probe 
of spin dynamics and equilibration. We first measure 
linewidth 7 of the mj = +1 transition in the dispersive 
regime [8] , with ~ 70 MHz ^> g ens detuning between the 
A/4 mode and mj = +1 transition. We extract 7 by 
inferring [23] the frequency shift (A/) of the A/4 mode 
immediately following a pump pulse whose frequency w p 
is stepped through resonance with the mj = +1 tran- 
sition [Fig. 4(b)]. The pump pulse slightly decreases 
the polarization of the ensemble, red-shifting the res- 
onator. We fit a Lorentzian lineshape to |A/|, finding 
a full-width-at-half-maximum 7/2tt = 9.0 ± 0.3 MHz. A 
similar dispersive measurement using the 3A/4 mode at 
B\\ = 263 mT gives 7/2tt = 12.0±0.7 MHz. We find these 
values to be temperature independent in the range 0.25- 
1.2 K [Fig. 4(c)], indicating that 7 is limited by dipolar 
interactions and field inhomogeneity. We attribute the 7 
increase with B» to the latter. 

The spin relaxation time is measured by applying 
a pump pulse resonant with the mj = +1 transition 
and monitoring the frequency shift in time as the spin 
polarization returns to equilibrium. We observe a bi- 
exponential decay response with time constants ~ 20 s 
and ~ 160 s [Fig. 4(d)]. These constants are independent 
of temperature in the range 0.25-1 K [Fig. 4(e)], suggest- 
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FIG. 4: (color online) Measurement of the spin linewidth 
and relaxation times using dispersive spin-resonator interac- 
tions, (a) Scheme (not to scale) for (b) the measurement of 
the spin linewidth (T = 0.25 K, By = 86 mT) by probing 
the frequency shift of the A/4 mode after applying a pump 
pulse (0.4 s duration, —50 dBm incident power) through res- 
onance with the mi = +1 transition (60 s wait between suc- 
cessive measurements) [23]. A similar measurement of 7 at 
By = 263 mT is obtained using the 3 A/4 mode, (c) 7 at 
B11 = 86 mT (circles) and Bu — 263 mT (squares) as a func- 
tion of temperature, (d) Measurement of the spin relaxation 
time Ti by probing the resonator shift as a function of time 
after the pump pulse is switched off. A bi-exponential decay 
is observed, (e) Temperature dependence of the two time con- 
stants, extracted by probing with the A/4 (circles) and 3A/4 
(squares) modes. Error bars, unless shown, are smaller than 
the symbol size. 



ing that spin polarization decay is not limited by spin- 
lattice relaxation [21] but spin diffusion instead. Through 
dipolar flip-flop processes, the depolarization diffuses out 
of the resonator mode volume, leading to repolarization 
of the ensemble. The rate for this process depends on the 
nominal dipolar coupling strength between spins in the 
ensemble, which itself depends on the spin density [12]. 
The two time constants may be explained by two dia- 
mond sectors inside the mode volume with electron-spin 
densities differing by a factor of ~ 8 [23, 26, 27]. This is 
supported by the fact that relative amplitudes of the two 
exponentials do not change versus temperature. 

To investigate spin dynamics across sub-ensembles, 
we measure how pumping one sub-ensemble can affect 
the coupling strength of other sub-ensembles to the 
resonator [23]. As shown in Fig. 5(a), pumping at 
/ mj =o(-Bii) completely suppresses the avoided crossing 
between the mj = transition and the resonator [28]. 
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FIG. 5: (color online), (a) Transmission spectroscopy similar 
to Fig. 2(b), with an additional pump pulse resonant with 
the mi = transition (incident pump power P p — —50 dBm, 
100 ms duration) prior to IS21I measurement. A complete 
disappearance of the mi = avoided crossing and a reduc- 
tion in the coupling strength of the undriven transitions are 
observed. Color scale is identical to Fig. 2(b). (b) Inset: 
vacuum-Rabi-split dips at B« = 269.1 mT as a function of 
Pp. The merging of dips with increasing P p indicates cross- 
relaxation between the sub-ensembles. Main panel: Extracted 
polarization P (normalized to the value without pump) for 
the undriven mi — + 1 sub-ensemble. The curve corresponds 
to the steady-state solution of a rate equation modeling fast 
equilibration between the sub-ensembles compared to Ti (see 
text and [23] for details). 



and polarizing magnetic field. The observed tempera- 
ture independence of spin linewidth and relaxation in the 
range 0.25-1 K supports spin out-diffusion as the dom- 
inant relaxation mechanism within the resonator mode 
volume. Resonant pumping of spin sub-ensembles indi- 
cates exchange of Zeeman and dipolar energies between 
sub-ensembles [20]. Follow-up experiments will pursue 
two directions: probing sub-ensemble response to one or 
more resonant pump pulses on millisecond timescales to 
shed light on the cross-relaxation mechanism, and cool- 
ing to 15 mK to fully polarize the ensemble [17] and 
extend spin coherence using sub-ensemble-selective dy- 
namical decoupling [18]. Ultimately, cross-relaxation and 
achieved coherence will set the timescale over which sub- 
ensembles may serve as independent quantum memories. 
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for Fundamental Research on Matter (FOM) , the Nether- 
lands Organization for Scientific Research (NWO, VIDI 
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Remarkably, partial depolarization is evident in the mj = 
±1 sub-ensembles. The coupling strengths of the un- 
driven transitions (mj = ±1) to the 3A/4 mode are re- 
duced to g cns /2TT = 12.5±0.5 and 12.0±0.5 MHz, respec- 
tively. To quantify this steady-state cross-relaxation, we 
measure the minimum-splitting between the hybridized 
dips at -By = 269.1 mT [arrow in Fig. 5(a)] as a function 
of pump power P p . As shown in the inset of Fig. 5(b), 
the undriven mi = +1 sub-ensemble depolarizes further 
with increasing P p . We can reproduce [23] this power- 
dependent steady-state cross-depolarization using a rate 
equation including a spin diffusion rate r o across the 
mode volume and a cross-relaxation rate T between sub- 
ensembles [20]. We assume T ^> T consistent with pre- 
vious measurements of cross-relaxation in high density 
Pl-center samples by Sorokin et al. [27]. Under these 
assumptions, the steady-state normalized polarization of 
each sub-ensemble is P = r o /(r o + O /3), where O is 
the pumping rate for the mi = transition. Excellent 
agreement is found with the model, with only the lever 
arm between Qq and P p as free parameter. Using the 
best-fit lever arm in combination with Fermi's golden rule 
£1 = 2ng 2 N phot /j and the measured r o w 0.05 s _1 and 
j/2tt 12 MHz, we estimate g ~ 2.5 Hz [29]. Compar- 
ing this g to g ens (T = 0) suggests N ~ 10 14 spins in the 
resonator mode volume. 

In conclusion, we have used resonant and dispersive 
interactions with the two lowest-frequency modes of a 
NbTiN CPW resonator to probe the dynamics of a PI 
electron-spin ensemble in diamond at low temperature 



[1] J. H. Wesenberg, A. Ardavan, G. A. D. Briggs, J. J. L. 

Morton, R. J. Schoelkopf, D. I. Schuster, and K. Mplmer, 

Phys. Rev. Lett., 103, 070502 (2009). 
[2] D. Marcos, M. Wubs, J. M. Taylor, R. Aguado, M. D. 

Lukin, and A. S. S0rensen, Phys. Rev. Lett., 105, 210501 

(2010). 

[3] X. Zhu, et al, Nature, 478, 221 (2011). 

[4] Y. Kubo, et al, Phys. Rev. Lett., 107, 220501 (2011). 

[5] H. Malissa, D. I. Schuster, A. Tyryshkin, A. Houck, , and 

S. Lyon, arXiv:1202.6305 (2012). 
[6] Y. Kubo, et al, Phys. Rev. B, 86, 064514 (2012). 
[7] A. Imamoglu, Phys. Rev. Lett., 102, 083602 (2009). 
[8] R. Amsiiss, et al, Phys. Rev. Lett., 107, 060502 (2011). 
[9] Y. Kubo, et al, Phys. Rev. Lett., 105, 140502 (2010). 
[10] D. I. Schuster, et al, Phys. Rev. Lett., 105, 140501 
(2010). 

[11] J. H. N. Loubser and J. A. van Wyk, Rep. Prog. Phys., 

41, 1201 (1978). 
[12] J. A. van Wyk, E. C. Reynhardt, G. L. High, and I. Ki- 

flawik, J. Phys. D: Appl. Phys., 30, 1790 (1997). 
[13] M. V. G. Dutt, L. Childress, L. Jiang, E. Togan, J. Maze, 

F. Jelezko, A. S. Zibrov, P. R. Hemmer, and M. D. Lukin, 

Science, 316, 1312 (2007). 
[14] P. Bushev, A. K. Feofanov, H. Rotzinger, I. Protopopov, 

J. H. Cole, C. M. Wilson, G. Fischer, A. Lukashenko, 

and A. V. Ustinov, Phys. Rev. B, 84, 060501 (2011). 
[15] M. U. Staudt, et al, J. Phys. B: At. Mol. Opt. Phys., 45, 

124019 (2012). 

[16] J. Clarke and F. K. Wilhelm, Nature, 453, 1031 (2008). 
[17] S. Takahashi, R. Hanson, J. van Tol, M. S. Sherwin, and 

D. D. Awschalom, Phys. Rev. Lett., 101, 047601 (2008). 
[18] G. de Lange, T. van der Sar, M. Blok, Z.-H. Wang, V. Do- 

brovitski, and R. Hanson, Sci. Rep., 2, 382 (2012). 
[19] R. Barends, et al, Appl. Phys. Lett., 97, 023508 (2010). 



5 



[20] N. Bloembergen, S. Shapiro, I. S. Pershan, and J. O. 

Artman, Phys. Rev., 114, 445 (1959). 
[21] E. C. Reynhardt, G. L. High, and J. A. van Wyk, J. 

Chem. Phys., 109, 8471 (1998). 
[22] R. J. Cook and D. H. Whiffon, Proc. R. Soc. Lond. A, 

295, 99 (1966). 
[23] See supplemental material. 

[24] E. Abe, H. Wu, A. Ardavan, and J. J. L. Morton, Appl. 

Phys. Lett., 98, 251108 (2011). 
[25] K. Sandner, H. Ritsch, R. Amsiiss, C. Roller, T. Nobauer, 

S. Putz, J. Schmiedmayer, and J. Majer, Phys. Rev. A, 

85, 053806 (2012). 
[26] R. C. Burns, V. Cvetkovic, C. N. Dodge, D. J. F. Evans, 



M. L. T. Rooney, P. M. Spear, and C. M. Welbourn, J. 
Crystal Growth, 104, 257 (1990). 

[27] P. P. Sorokin, G. J. Lasher, and I. L. Gelles, Phys. Rev., 
118, 939 (1960). 

[28] This is surprising from a single-spin perspective, be- 
cause the maximum Rabi driving strength (/fum = 
g^/N phot /2iv « 100 KHz for Ar phot = 10 8 photons on 
mode resonance) is significantly smaller than the spin 
linewidth. 

[29] Note that iVphot is lower than on mode resonance by the 
filter factor (k; + k c ) 2 /(cj c — oj p ) 2 . 



Supplement to "Probing dynamics of an electron-spin ensemble via a superconducting 

resonator" 



V. Ranjan*, 1 G. de Lange*, 1 R. Schutjens, 1 T. Debelhoir, 2 J. P. Groen, 1 

D. Szombati, 1 D. J. Thoen, 1 T. M. Klapwijk, 1 R. Hanson, 1 and L. DiCarlo 1 

1 Kavli Institute of Nanoscience, Delft University of Technology, 
P.O. Box 5046, 2600 GA Delft, The Netherlands 
2 ICFP, Departement de Physique de I'ENS, 24 rue Lhomond, 75005 Paris, France 

(Dated: August 28, 2012) 



MEASUREMENT SETUP 

The device is operated in a 3 He refrigerator with base 
temperature T = 0.25 K. All measurements involve 
continuous-wave heterodyne detection of feedline trans- 
mission Sol function of Bu and T with 1 MHz in- 
termediate frequency and 80 ms integration. Incident 
power on the feedline ranges from —115 to —95 dBm 
(-110 to -90 dBm) for the A/4 (3A/4) mode, corre- 
sponding to ~ 10 3 to 10 5 intra-cavity photons on res- 
onance (incident power for one photon on resonance 
w hjj c (ni + k c ) 2 /2k c [1]). All other microwave pulses 
are also applied through the feedline. At B\\ — the 
diamond-coupled resonator has A/4 and 3 A/4 modes at 
2.59 and 7.67 GHz, respectively (~ 11% lower than the 
values without the diamond). The reduction in fre- 
quency of the resonator after mounting diamond is due 
to the higher effective dielectric constant of the diamond 
(<~ 5.7). The corresponding loaded quality factors are 
Q = w c /(/Ci + K e ) = 1.8 x 10 4 and 6 x 10 3 (~ 3 times 
lower than the values without the diamond). Increase in 
intrinsic losses could come from the surface of the dia- 
mond and from the vacuum grease used to glue the dia- 
mond to the resonator. Increasing Bn to 300 mT further 
decreases Q by ~ 15%. To reduce microwave losses in 
the superconducting film resulting from trapped flux, the 
system is frequently reset by warming above the film T c 
and cooling back down with B\\ = 0. 

AVOIDED CROSSINGS 

Transmission spectroscopy measurements in Fig. 2 re- 
veal two regimes of spin-resonator coupling. Rabi-split 
dips are clearly resolved [Fig. 3(a)] when the mj = tran- 
sition is resonant with the 3 A/4 mode, evidencing strong 
coupling (2g ons > 7, k). To extract g cns in this case, we fit 
Eq. (1) to IS21K/) [2]. For this fit, we fix m, k c and 7 to 
independently-measured values: k;/27t = 1.4 ± 0.1 MHz 
and k c /2tt — 120 ± 10 KHz are extracted on both sides 
of the avoided crossings, and 7/27r = 12.0 ± 0.7 MHz is 
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obtained from dispersive pump-probe linewidth measure- 
ments (Fig. 4). 

The A/4 mode is weakly coupled to the resonator 
(2<?ens < 7) and Rabi-split dips are not resolved on 
resonance. The extraction of a coupling strength to 
this mode involves several steps. First, we extract the 
field-dependent loaded quality factor Q by fitting a line- 
shape [1] 

S£ a + i2Q(u-Lj I )/u, I 
21[UJ) 1 + i2Q(uj-u i )/uj i 

to the measured jS^iKw) at each magnetic field. Fig- 
ure 3(b) shows Q(B\\) near the avoided crossings for the 
A/4 mode at T = 0.25 K. The decrease in Q when 
the spin transitions are tuned into resonance with the 
resonator is caused by absorption of microwave energy 
by the spin ensemble. An analytic approximation of 
Q = o; r /Acj r valid away from the avoided crossings can 
be obtained from Eq. (1) and noting that on resonance, 
where |52i| is minimum, Im(S2i) ~ 0. Here, 

A^ r = («i + K C ) + flens (A 2 ^ 2/4) 

and A n = cj r — uj mi=n . The fixed values, 7/2-7T = 9.0 ± 
0.3 MHz, Ki/27r = 140±10 KHz and k c /2it = 32±1 KHz, 
are obtained by the same methods as described for the 
3A/4 mode. 

The extraction of g C ns requires accounting for all en- 
sembles in the last term of Acj r . We also observe two 
small dips ~ 4.1 mT away from the avoided crossings 
with the mi = ±1 transitions [indicated by arrows in 
Fig. 3(b)]. These dips are consistent with coupling to the 
nuclear spin (I = 1/2) of a 13 C (1.1% natural abundance) 
at the nearest-neighbor carbon site present on the Jahn- 
Teller axis (Ai3 C ^ = diag(141.8, 141.8, 340.8) MHz [3]). 
These give rise to six sub-ensembles with coupling 
strengths g^ ns i 3c ~ ffens ^0.011/2. The factor 1/2 comes 
from the fact that 13 C has nuclear spin of I = 1/2 in con- 
trast to 14 N (7 = 1) Additionally, there are six more C 
transitions with Ai3 Cib = diag(32.1, 32.1, 41.0) MHz [3], 
which arise from 13 C located at the three remaining 
nearest-neighbor sites. These are three times more abun- 
dant with Sen S ,i3C ~ V^Sens^C- Although they are 
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not visible as separate dips in the data, they do modify 
Q(Bu). We therefore include all of these sub-ensemblcs 
in our modeling. 

One can in principle also find an analytical expression 
for the dispersive shift of the resonance with respect to 
ui c and use it to extract g ens for the A/4 mode. In our 
experiments, however, this is problematic due to sudden 
changes in cj c [indicated by the arrow in Fig. 2(a)] caused 
by flux jumps in the NbTiN film. The absorption, how- 
ever, only depends on w r , which is extracted with high 
precision from the transmission lineshape. 



DISPERSIVE MEASUREMENTS 

When performing the dispersive measurements shown 
in Fig. 4, we measure the change in IS21I at fixed fre- 
quency / = (w r — Aw r /2) /2n. Immediately before the 
pumping step, we measure the resonator lineshape and 
calibrate the slope of IS21I at /. We use this calibration 
to infer the resonator shift A/. To ensure linearity within 
10%, the pump pulse power and duration are chosen as 
to keep the maximum |A/| < Aw r /87r. 

Measuring the relaxation of pumped spins over time re- 
veals bi-exponential decay curves. A possible reason for 
this observation would be the presence of two sectors with 
different spin densities within the mode volume. This is 
suggested by several tests. First, we repeat the exper- 
iment with the pump far detuned from any spin tran- 
sition and observe no shift in the resonance frequency. 
This excludes the possibility that one of the time con- 
stants results from the pump influencing the resonator 
shift (for example, through heating). Second, we per- 
form the experiment varying the pump pulse frequency 
around f mi=+ i and observe that the amplitudes of both 
exponentials are modulated by the spin linewidth (see 
Fig. SI). This excludes the possibility of a broad back- 
ground spin ensemble such as paramagnetic impurities on 
the surface of NbTiN or diamond itself. In addition, the 
ratio of g cns (T — 0) values for the two modes differs from 
the expected V% by a factor ~ 1.6. This further hints 
at a strong inhomogeneity of PI concentration across the 
diamond. 



PUMPING AND DEPOLARIZATION OF 
ENSEMBLES 

To study the dynamics of spins within and between 
the sub-ensembles, we resonantly pump the mj = 
transition and measure the effect on the sub-ensembles. 
The pulse sequence used for transmission spectroscopy 
|S2i|(/,-B||) in Fig. 5(a) is shown in Fig. S2. 

We model the depolarization of spin ensembles due to 
pumping with a simple rate equation. In addition to 
the pumping rate £Iq for the driven transition (to/ = 0), 
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FIG. SI: Best-fit amplitudes of the two exponential terms 
observed in the relaxation measurement as function of pump 
frequency. Lorentzian fits have widths "fi/2n = 8.5±0.5 MHz 
(solid circles) and 72/2^ = 10.7 ± 1.5 MHz (open circles) for 
fast and slow decays, respectively. Within error bars, these 
widths are consistent with the spin linewidth measurements 
in Fig. 4(b-c). 



we include a spin-equilibration rate r o (spin diffusion 
rate across the resonator mode volume) within each 
sub-ensemble and a nominal rate T that characterizes 
cross-relaxation between the sub-ensembles. The rate of 
change of polarization (equal for mj — ±1 due to sym- 
metry), is modeled by 



d_ 
dt 



with 



P = 



-MP + r, 



p(±iy 

P(0) 



(SI) 



and 



M = 



r = r G 



r Q + r 
-2T r n 



-r 

O - 



2r 



The steady-state polarization (dP/dt = condition) is 
P = M~ 1 T, yielding 



P(±l) 



r o (r o + 3r + » ) 
r o (r o + 3r + fi ) + r^ ' 



For r r o [4], this expression simplifies to P(±l) = 
P(0) = r o /(r o + O /3). In this limit, the fast equili- 
bration between sub-ensembles causes pumping of each 
transition with an effective rate £7/3, even though only 
one transition is directly driven by the pump. 
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FIG. S2: Schematic of pump- and probe-tone timings used in 
Fig. 5. At each B\\ setting, following a 5 s waiting time, pump 
pulses resonant on driven transitions are applied for 20 s. This 
is followed by a periodic series of pumping and probing steps 
(period 2r), with the probe frequency / stepped at the start 
of every cycle. The probe is always on, but l&iK/) only 
measured for a time r — r'. Here r' is a settling time of 
20 ms. We use r = 100 ms for 2-D transmission spectroscopy 
[Fig. 5(a)], and r = 270 ms for measurements of Rabi-split dip 
separation as a function of pump power at fixed B« [Fig. 5(b)]. 



